Band-filling effects on the Kondo-lattice properties 
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We present theoretical results for a Kondo-lattice model with spin- 1/2 localized moments, in- 
cluding both the intrasite Kondo coupling and an intersite antiferromagnetic exchange interaction, 
treated within an extended mean-field approximation. We describe here the case of a non-integer 
conduction-band filling for which an "exhaustion" problem arises when the number of conduction 
electrons is not large enough to screen all the lattice spins. This is best seen in the computed mag- 
netic susceptibility. The Kondo temperature so obtained is different from the single-impurity one, 
and increases for small values of the intersite interaction, but the Kondo-effect disappears abruptly 
for low band filling and/or strong intersite coupling; a phase diagram is presented as a function of 
both parameters. A discussion of experimental results on cerium Kondo compounds is also given. 

PACS numbers; 71.27.-|-a, 75.30.Mb, 75.20.Hr, 75.10.-b 



I. INTRODUCTION 

The properties of many cerium or ytterbium com- 
pounds are well accounted for by the Kondo-lattice 
model, where a strong competition exists between 
the Kondo effect and magnetic ordering arising from 
the RKKY (Ruderman-Kittel-Kasuya-Yosida) interac- 
tion between rare-earth atoms at different lattice sites- 
This situation is well described by the Doniach diagram,Ej 
which gives the variation of the Neel temperature and 
of the Kondo temperature with increasing antiferromag- 
netic intrasite exchange interaction Jk between local- 
ized spins and conduction-electron spins. If one consid- 
ers the exchange Hamiltonian between localized (S) and 
conduction-electron (s) spins, given by 



H = Jks-S, 



(1) 



usual theories of the one- impurity Kondo effect and of the 
RKKY interaction yield a Kondo temperature Tko that 
is proportional to exp(— l//9Jx), and an ordering temper- 
ature (Neel or, in some cases. Curie) Tno, proportional 
to pJ^, p being the density of states for the conduction 
band at the Fermi energy. Thus, for small pj^ values, 
Tno is larger than Tko and the system tends to order mag- 
netically, with often a reduction of the magnetic moment 
due to the Kondo effect. On the contrary, for large pJk, 
Tko is larger than Tno and the system tends to become 
non magnetic. The actual ordering temperature T/v, 
therefore, increases initially with increasing pJk, then 
passes through a maximum and tends to zero at a critical 
value pJ^ corresponding to a "quantum critical point" 
(QCP) in the Doniach diagram. Such a behavior of Tn 
has been experimentally observed with increasing pres- 
sure in manyijerium compounds, such as CeAl2 (Ref. |^) 
or CeRh2Si2.a We also know that the Neel tempera- 
ture starts from zero at a given pressure, and increases 
rapidly with pressure in YbCu2Si2 (Ref. ^ or in related 
ytterbium compounds, which can be considered as an- 
other test of the Doniach diagram. The one-impurity 



model predicts an exponential increase of the Kondo 
temperature with pJk- This means that the Kondo 
temperature should increase with increasing pressure in 
cerium compounds and with decreasing pressure in ytter- 
bium compounds, in good agreement with many observa- 
tions. However, deviations seem to occur in some cerii. 
compounds, such as CeRh2Si2 (Ref. ^ or CeRu2Ge2j 
where the actual Kondo temperature observed in a lat- 
tice can be significantly different from the one derived for 
the single-impurity case. Thus, in order to account for 
such an effect, we have treated in a previous paperu the 
Kondo-lattice model with both intrasite Kondo exchange 
and intersite antiferromagnetic exchange interactions, for 
a half-filled conduction band (corresponding to a number 
of conduction electrons n — 1). We employed a mean- 
field approximation with two correlators, Aj^, describing 
the intrasite Kondo correlation, and r^^-^, representing an 
intersite correlation between two neighboring moments. 
We have shown that the enhancement of the intersite ex- 
change interaction tends to decrease the Kondo temper- 
ature Tk for the lattice with respect to the one-impurity 
Kondo temperature Tko, and to suppress the Kondo ef- 
fect for large values of the intersite exchange interaction 
parameter.13 So, this model can account for the pressurje 
dependence of Tk observed in CeRh2Si2j3 CeRu2 Ge2 Jj'EI 
or more recentlyB Ce2Rh3Ge5. 

The above mentioned model also yields a correla- 
tion temperature Tcor, below which short-range mag- 
netic correlations between neighboring cerium atoms oc- 
cur. For sufficiently large values of the intersite ex- 
change, Tcor is higher than the Kondo temperature. This 
result is in agreement with the experimental observa- 
tion by neutron diffraction experiments of such shorty 
range magnetic correlations in single crystals of CeCu6,u 
CeInCu2,t3 CeRu2Si2MEiy or Cei_^La^Ru2Si2 (Refs. 
|ll|,|l2 13 1^) at low temperatures. The experimentally 
observed temperature Tcor is clearly larger than the 
Kondo tempofl^^e T^: Tcor-- 60-70 K and Tk~ 20 



in CeRu2Si2 



Tcor-- 10 K and T^- 5 K in CeCue.l 
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Most Kondo-lattice models have been studied for the 
case of a half-filled conduction band, corresponding to 
a number of conduction electrons, n, equal to the num- 
ber Hf — 1 of f electrons. However, when n < 1, an 
"exhaustion" problem arises, which means that there are 
not enough conduction electrons to screen all the local- 
ized gpins and as a consequcncc, the Kondo temperature 
decreases. The stability of the Kondo effect when n 
decreases has been recently studied,tZl and also the tem- 
perature dependence of the magnetic j^sceptibility has 
been computed for different values of ri,Ea with the Hamil- 
tonian including an intersite exchange interaction Jh be- 
tween the localized spins, treated as independent of Jk- 
The results showed that a reduction of n and/or an en- 
hancement of Jh tend to suppress the Kondo effect. In 
particular, the analytical calculation at_^ — gave 
proportional to n in the case of Jh = 0£3 

Early studies of the Kondo-lattice model h a v e , a on- 
sidered only the intrasite Kondo interaction.E30n3 In 
particular, Continentino et al.jHd using scaling theory, 
have found a coherence temperatuxp.. increasing above 
the QCP. Recently, Burdin et al.Oo using the slave- 
boson method in the exhaustion limit, have obtained 
two different energy scales: the single-impurity Kondo 
temperature Tk, which corresponds to the onset of lo- 
cal singlet formation, and the zero-temperature energy 
gain, T*, related to the coherent Kondo effect. Thus, the 
low-temperature energy scale is not equal to the high- 
temperature one, altho ugh t here is only one mean-field 
parameter in the thcory.E^O It has been shown that the 
ratio T* /Tk is much smaller than 1, and depends only 
on the band filling factor. In the case of a rectangular 
density of states for the conduction band, Tk vanishes 
as v}/"^ and the coherence temperature as n in the limit 
n ^ 0. However, both Tk and T* where found to ex- 
hibit the same exponential behavior exp(— 2D/J«), in 
contrast with the formula established by Nozieresjl3 who 
found T* ~ T^/ D. The effect of a small number of con- 
duction electrons has been also studied within. both the. 
Kondo-lattice and the Anderson-lattice models.L3t3E!l'E3 

The purpose of the present paper is to further de- 
velop the study of conduction-band filling effects within 
the KandcJattice model. Besides previous work cited 
aboveJ^ 



a first brief account has been already 
given,Ea but we will present here both numerical results 
and new analytical expressions for some specific limits, 
as well as a detailed discussion of band-filling effects. 



II. THE THEORETICAL MODEL 

As it was mentioned above, the Kondo-lattice model 
is often studied with only the intrasite Kondo interac- 
tion. In principle, both the RKKY magnetic interaction 
and the Kondo effect are consequences of the intrasite 
exchange term. However, the RKKY interacticm, is per- 
turbative in Jk while the Kondo effect is not.E3 When 
dealing with approximation schemes, it is hard to obtain 



both effects starting with the intrasite interaction alone. 
For this reason, an explicit intersite exchange intctaction 
is usually included in the model Hamiltonian.Clca'L3 
The proposed Hamiltonian of the system is, therefore. 



H 



Jk ^ Si • Si -I- Jh ^ • 



(2) 



where Sk is the energy of the conduction band, Jk is 
the Kondo coupling between a localized spin Si and the 
spin Si of a conduction electron at the same site, and Jh 
is the interaction between nearest-neighboring localized 
spins. Assuming spin-1/2 localized moments, we repre- 
sented them by a zero-width / band with one electron 
per site, while the conduction band has width 2D and a 
constant density of states. We choose Jk and Jh to be 
positive, implying that both local and intersite interac- 
tions are antiferromagnetic, as is the case in most cerium 
compounds. 

We now write the spin operators in fermionic represen- 
tation: 



t = c^Va , St = fy.i , (3) 



remembering that we have a constraint of single- 
occupancy of the / level at all sites, n{ = 1. 

In order to discuss the Kondo effect and magnetic cor- 
relations we define the operators 



A,-, 



f if 



V f 

J 1(7 J J C 



(4) 



where Ai^ describes the intrasite Kondo correlation, 
and ^ijcr represents an intersite correlation between two 
neighboring atoms. This is so because the on-site spin- 
spin correlation (si • Si) can be written in terms of the 
average {\^^\_^), and similarly for the intersite cor- 
relation using f j^g.. With this notation we perform an 
extended mean-fiikl approximation, introduced by Cole- 
man and Andrei,E2l and presented in full detail in Ref. |^. 
Considering translational invariance, and taking into ac- 
count that there is no breakdown of spin symmetry, i.e., 
no magnetic states, we can write A = (Aio-) for all sites, 
and r = (Xijiy) f'^r nearest-neighboring sites and zero oth- 
erwise. In this way we obtain a mean-field Hamiltonian 
that takes the form of a hybridized two-band system: 

kcr i \ a / 

- J^,A^(4/,,,-f /il^cj 



with 



E (/-/j- + 4/-) + , (5) 



Ek = 2NJkX^ , Eh = zNJhT^ , (6) 
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TV being the total number of lattice sites. We have intro- 
duced a term depending on the Lagrange multiplier Eq 
in order to impose the constraint '^i{n( — 1) = 0. This 



/ _ 



1. 



is a weak form of the original constraint n 

After performing this approximation, one deals with 
a one-electron Hamiltonian representing two hybridized 
bands: the conduction band of width 2D and the / band 
of effective band width 2BD, with B given by 



B = -zJhT/D, 



(7) 



z being the number of nearest neighbors of a site, while 
the magnitude of the hybridization gap is directly related 
to . This quantity is also a measure of the Kondo effect, 
as the Kondo correlation function (s^ • S^) is proportional 
to A^. 

The Hamiltonian (|^), leaving aside the constant terms, 
is easily diagonalized, and the resulting energies of the 
two new hybrid bands are 



ei,{l + B)+Eo 



4Ji.A2 



(8) 



The mean-field parameters A and T are obtained by 
self-consistently solving Eqs. (^) or, equivalently, by min- 
imizing the total internal energy 

E^2 E^f{E^)+EK + EH^EoN (9) 

k,Q=± 

at zero temperature, or the total free energy F at finite 
temperatures .Il3 The summation in Eq. is made over 
all k states in the first Brillouia zone, and /(i?k) is the 
Fermi-Dirac function. As usual,l3 the reference energy Eq 
of the / band and the chemical potential ^ have to be 
determined self-consistently in order to keep the average 
numbers of / and conduction electrons respectively equal 
to 1 and n. 

We have previously treated the half-filled case, corre^ 
sponding to a number of conduction electrons n = 1,0 
and here we consider the general case of a non-intege* 
number of conduction electrons. For the half-filled case,Ll 
we have shown that the Kondo temperature Tk for the 
lattice can be much smaller than the Kondo temperature 
Tko for the single impurity, as observed in some cerium 
compounds. We will present in the next sections the the- 
oretical results for band-filling effects, and discuss later 
on the comparison with experiment. 



III. BAND-FILLING EFFECTS AT T = 



The two hybridized bands E^ given by Eq. (| 
a structure that depends on the two factors A 
and B, defined by (|^), and especially on the sign of B. In 



exhibit 




FIG. 1: Schematic plot of the non-hybridized conduction and 
4/ bands, and of the two hybridized bands E± (k) . For a given 
n < 1, the chemical potential is different in panels a) and 
b), which correspond respectively to the two cases studied; 
small Jh (a single solution of _E_ {e) = ^) , and large Jh (two 
solutions). All energies appear in units of D. 



the case of n = 1, previously considered ,0 for small \B\ 
values such that |i?|_D^ < A there is a gap, the lower band 
is completely full and the upper band empty at T = 0. 
For n < 1, the Fermi level cuts the lower band, and the 
upper band is still empty at T = 0, which corresponds, 
therefore, to a real metallic situation. 

The shape of the lower band E^^ and the solution of 
the case n < 1 depend critically and self-consistently on 
the values of the different parameters Jk, Jh and n or 
equivalently A, B and n. li B > 0, i.e., F < (since Jh 
is positive here), both the "/ band" and the hybridized 
lower band E^ are continuously increasing with e^. The 
same structure exists for i3 < and small Jh values, 
while E^ presents a maximum for B < and larger Jh 
values. 

Thus, the shape of E^ can change under a variation 
of the parameters, and this peculiar situation makes the 
problem difficult to solve. The two characteristic cases 
are presented schematically in Fig. |l|, where we have used 
the fact that the energies E^ depend on k only through 
the bare conduction-band energies £(k), making the sub- 
stitution E^ —> E±{e), with e defined in the interval 
[-1,1] in units of the half bandwidth D. As can be seen 
in Fig. |l|, the curve (e) may or may not have a maxi- 
mum, according to the different cases studied. Moreover, 
when E- (e) presents a maximum, the Fermi energy can 
cut the lower band in two points or in only one point 
within the first Brillouin zone. In fact, we have to distin- 
guish two cases, depending on whether there are one or 
two intersection points (which implies automatically in 
the second case that the curve £'-(e) has a maximum). 
We will compute the total energy in both situations. 

For T = 0, we firstly derive analytical expressions 
for both cases, and then we present numerical solutions 
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for different sets of parameters that affect the shape of 
E_{e). In the first case, shown m Fig. Q a), the energy 
per lattice site is written as 



give the final values of A and T corresponding to the 
minimum energy. For small values of Jh , Eq. ( |l3| ) allows 
us to write 



^ r de 
2D J_o 



Eo + il + B)e 



(10) 



The upper limit Eq is the value of Ek that corresponds 
to E^ = ^ in Eq. (|8|), as shown in Fig. |l|. Writing the 
two self-consistent equations that yield n/ = 1 and the 
number of conduction electrons n, we obtain — Dn, 
and 



with 



Eo^--il-n) [il-B)D + A] 



(11) 



(12) 



Finally, the internal energy per site at T = [Eq. (10)] 
is given by 



J|.A2 f A - 7:>(1 - B) 



D{l-B) \/\ + D{1- B) 



(13) 



The derivatives of the total energy with respect to A 
and F yield two self-consistent equations, whose solutions 



A^ 
F 



nD'^il - B) 

nA 
4i:>(l - B) 



[l-uf 
l-n 



D{l-B) 



(14) 
(15) 



with u = exp [-2D{1 - B)/Jk]- 

The case with a maximum in E^ is more complicated. 
As shown in Fig. |l|, for the case B < with large values 
of Jh, the Fermi level fj, cuts the band E^ in two ener- 
gies; then, the two corresponding energies ei and S2 are 
given by the values of e(k) which are roots of the equa- 
tion E^^ — /i with E^ given by eq.(H). So, the expression 
of the energy is given by Eq. ) by changing the inte- 
gration limits so that the Integra is now performed from 
~D to ei and from £2 to D. 

Now, the two equations derived from the constraint 
Uf — 1 and the number of conduction electrons equal to 
n are 



£2 - £1 = D{1 - n) , 



Eo = -{l-n)JD'^ 



4JiA2 



(l-B)2-(l-n)2 



(16) 
(17) 



We can derive analytical expressions for £1 and £2, and 
finally for the total energy per site £ for _B < 0. A quite 
long calculation yields 



J 



R 



Ail - Bf - 2B{1 - n) + (I - -n^)] 
2(1 - B) ' ^ ' ^ {l-B)D 



In. 



P - D{l-n)^/\B\ R + D 
P + D{l-n)^/\B\ 



2(1 -s; 



■P\/\B 



D'B' 
zJh 



2JkX^ — Eq 



(18) 



where 



R 



4JiA2 



(l-S)2-(l-n)2 



(19) 



and 



P = ^^4JlX^ - D^B{l-n)'^ . (20) 
Minimizing the energy (|l^) with respect to A we obtain 
A2=p/2a, (21) 



with 

P^2Jk 
and 



In 



D{l-B) 



(l-B)2-(l-n)2 



-B{l-ny 



(22) 



72 



AD{1 - B) 



B(l-7l)2 (l-B)2-(l-n) 



(23) 

For the present case, A vanishes at a critical band filling 
given by 



= 1 - (1 - B) 



u- B 



(24) 
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0.0 0.2 0.4 



FIG. 2: Variation of A and F with the band filhng n for 
Jk/D = 0.4 and different values of Jh- 



We have finally performed a numerical calculation by 
considering the general form of the energy given by either 
Eq. ( p^ or Eq. (|l^) in the different cases studied, and 
looking for the minimum energy of the system. In fact, 
for each value of the different parameters, we compare 
the two solutions and we retain the lowest-energy one, 
thus determining the corresponding A and T . 

Figure || shows and F versus n for different values 
of JH^ and Jk/D = 0.4. For Jh = we obtain negative 
values of F (not shown in Fig. |^), with F —1/4 in 
the limit n ^ 0- see Eq. (|l^)] . This is consistent with 
previous resultalII'E2l showing that F changes sign as a 
function of n when the Kondo effect is strongly dominant. 

Figure ^ shows the phase diagram for Jk /D — 0.4 and 
1.0. In the latter case (previously reported in Ref. 29), 
for small Jh the Kondo phase is stable for all values 
of the band filling, while for large Jh the Kondo phase 
is stable only for n > ric given by Eq. (p^). Figure ^ 
also shows that Eq. (|2^) is very close to the numerical 
result in this region, but does not apply when Jh is very 
small, since the assumed maximum in the lower energy 
band no longer exists. The crossover between these two 
regimes was obtained numerically. In the crossover region 
we found a discontinous transition from At^O to A = 0. 
These results have to be taken with caution, however, 
since the value of Jk/D = 1 is unphysically high, and the 
Kondo lattice is expected to show ferromagnetic behavior 
in the \aw-n limit J23 while our analysis is restricted to 
antiferromagnetic intersite exchange. 



Thus, our calculation shows clearly that small n and 
large Jh values tend to suppress the Kondo effect, yield- 
ing a "magnetic" phase with A = and large short-range 
magnetic correlations. In fact, in this region a long-range 
magnetic order should certainly be stabilized, but this 
was not taken into account in this approach. In contrast, 
both A and F are different from zero in the Kondo phase. 



FIG. 3: Phase diagram plotted as the critical band filling ric 
versus Jh- The curves are drawn from Eq. ( |24| ) for Jk/D = 
0.4 (dotted line) and 1.0 (solid line). The symbols correspond 
to the results obtained by minimizing the energy. In each case, 
the Kondo regime is stable above the line, and the magnetic 
phase below. 



IV. BAND-FILLING EFFECTS AT FINITE 
TEMPERATURE 

Focusing on the Kondo-lattice problem at finite tem- 
peratures, the half-filled case has been previously de- 
scribed, and we consider here the general case of n < 1. 
The number of / electrons n/ is always taken equal to 1. 
The values of A and F are determined by self-consistently 
solving Eqs. (^ or by minimizing the free energy, which 
is given by 



F = -2T ^ In 



1 + e 



Ea + zJhT"^ + 2Jk\^ 



(25) 



In our mean-field approximation, Tk and Tcor are de- 
fined as the temperatures at which respectively A and F 
become zero. We have shown, in the case n = 1, that 
the Kondo temperature Tk for the lattice can be much 
smaller than the Kondo temperature Tko for the single 
impurity, as observed in some cerium compounds J3 More- 
over, we have also shown that the temperature Tcor for 
the occurrence of short-range magnetic correlations lies 
above the Kondo temperature Tk for large Jh values, 
while Tcor = Tk for small Jh- These two theoretical re- 
sults have provided a good explanation for the previously 
described anomalous behavior observed in some cerium 
compounds. 

First, we present analytical results obtained in the case 
of large Jh values, where Tcor > Tk- In this case, Tcor 
is easily determined by taking the limit F — > in the 
self-consistent equation for F with fixed A = 0, and we 
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FIG. 4: Kondo temperature Tk versus Jh for Jk/D = 0.4 
and several values of n. We also show the correlation temper- 
ature Tcor- 



obtain 



zJh /12 , 



(26) 



where z is the number of nearest neighbors. In the simple 
cubic case, z = 6, we obtain Tcor = Jh/'^- Still for the 
case Tcor > Tk , we determine Tk by taking the limit A — > 
in the self-consistent equation for A, with T remaining 
finite. We obtain the following equation to determine 
Tk- 



1 



pJk 



de 



-D 



1 



e(l -B)-Eo 
1 



1 + 



1 



(27) 



Equation (^7|), together with the self-consistent evalua- 
tion of B [see Eqs. (|^) and (^)], and the two equations 
giving Uf = 1 and n < 1 (which yield Eq and fi) allow us 
to determine Tk- In the case of Jh ~ 0, we recover the 
simple equation 



pJk 



-D 



de 



-D 



1 



e-Eo 



■ tanh 



2Tk 



(28) 



which reproduces the single-impurity result. In the limit 
Jk/D « 1, this yields 



Tk = CqD^ n{2 — n) exp 



1 



pJk 



(29) 



where Cq is a numerical constant: Co = 1.1337. Eq. (E£ 



gives Tk depending on the band filling as \Jn(2 — n), 
which is consistent with a ^/rl dependence when n ^ 0, 
as found in Ref. ^ 

Now, we present numerical results obtained at finite 
temperatures in the case n < 1, for different sets of pa- 
rameters Jk, Jh and n, in order to study the Kondo- 
lattice properties and band-filling effects. Figure ^ gives 



FIG. 5: Kondo temperature Tk as a function of Jk for 
Jh /D — 0.04 and representative values of n. We also show the 
correlation temperature Tcor, and the single-impurity Kondo 
temperature Tko- 



the curves of the Kondo temperature Tk versus Jh for a 
given Jk value and several values of the conduction-band 
filling n. We see clearly that Tk first increases, and then 
decreases with Jh for fixed n, dropping abruptly to zero 
at some critical value of Jh- On the other hand, for a 
given Jh, Tk decreases rapidly as n departs from half 
filling. The enhancement of Tk for small Jh was masked 
in Ref. by the choice of Jh — aJ]^, which intended to 
mimic the RKKY interaction. Besides this being a ques- 
tionable approximation, we have preferred here to study 
the behavior of the Kondo temperature as a function of 
Jk and Jh considered as independent parameters. For 
comparison, in Fig. ^ we have also plotted the correla- 
tion temperature Tcor, which is linear with Jh indepen- 
dently of the value of n.Q Figure || gives Tk as a function 
of Jk for Jh/D = 0.04 and representative values of n. 
Here again we include the correlation temperature Tcor, 
which signals the onset of short-range magnetic correla- 
tions when the temperature is lowered at fixed Jk- For 
comparison, we plot the single-impurity Kondo temper- 
ature Tko (for n ^ 1), which varies exponentially with 
Jk- and is weakly dependent on n near half-filling [see 
Eq. (H)]. 

Figures ^ and ^ show some interesting results of our 
model. First, one can see the occurrence of short-range 
magnetic correlations above the Kondo temperature, in 
good agreement with experiment. Also, in the region of 
coexistence between Kondo effect and magnetic correla- 
tions, the Kondo temperature is enhanced with respect to 
the single-impurity case, and shows a smoother variatioa 
with Jk- Such an enhancement agrees with the ideaE2l 
that /-/ correlations may, in some cases, reinforce the 
Kondo effect on a lattice. The second noticeable feature 
of Figs. ^ and ^ is the almost catastrophic supression of 
the Kondo effect with increasing intersite coupling, and 
the enhancement of this behavior as the band-filling fac- 
tor is reduced. 



7 




0.025 



T/D T/D 



FIG. 6: Magnetic susceptibility times temperature as a func- 
tion of temperature for n = 0.8, Jk/D — 0.4, and different 
values of Jh- The open circles mark the position of Tcor, and 
the stars, that of Tk- These two temperatures coincide in the 
case of Jh/D = 0.02. 



FIG. 7: Magnetic susceptibility times temperature as a func- 
tion of temperature for Jk/D — 0.4, Jh/D — 0.04 and differ- 
ent values of n. The open circle signals Tcor, while the filled 
circles indicate the positions of Tk for the three higher values 
of n (the Kondo regime no longer occurs for n — 0.6). 



Now, in order to better understand the physics of 
the Kondo lattice at low temperatures, and to see more 
clearly the different regimes occurring with decreasing 
temperature, we present different curves showing the 
product of magnetic susceptibility and temperature, x T, 
which represents the square of the effective magnetic mo- 
ment. 

On the two following figures, x ^ is constant at high 
temperatures down to a first kink corresponding to ei- 
ther the correlation temperature Tcor for a sufficiently 
large Jh value or to Tk at low Jh values. The kink 
corresponding to Tk in the case when Tk < Tcor is less 
pronounced, but still clearly visible. Figure |^ gives the 
results for the temperature dependence of the product xT 
for Jk/D = 0.4, n = 0.8, and representative values of Jh. 
This figure shows a Kondo regime alone for small Jh, a 
regime with only short-range magnetic correlations above 
the Kondo regime for intermediate Jh, and the absence 
of the Kondo phase for large Jh- Then, Fig. |^ gives the 
results for the temperature dependence of the product 
xT for Jk/D = 0.4, an intermediate value Jh/D = 0.04, 
and representative values of n. The curves show the on- 
set of short-range correlations and the Kondo regime for 
n values not very far from 1, while the Kondo effect has 
completely disappeared for n — 0.6, in which case a sin- 
gle critical temperature remains. We immediately see 
on this figure that the behavior of the magnetic suscep- 
tibility is completely different for n — 1 and smaller n 
values, because the n = 1 case corresponds to an insulat- 
ing system, and smaller n values yield metallic behavior. 
We also see that the Kondo effect disappears when n de- 
creases below a critical value, which can be iriterpreted 
as a manifestation of the "exhaustion" problemllj in con- 
nection with the competition between Kondo effect and 
magnetic correlations. 

Finally, in figure 0, we present a plot of the magnetic 
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FIG. 8: Magnetic susceptibility as a function of temperature 
for Jk/D — 0.4 and different sets of the parameters Jh and 
n, in order to show the different regimes. Circles and stars 
follow the convention of Fig. H. 



susceptibility versus temperature for the three different 
regimes for typical sets of parameters: Kondo regime 
alone (small Jh, n near half- filling) , both magnetic cor- 
relations and Kondo regimes (higher Jh , smaller n) , and 
short-range magnetic correlations only (still higher Jh, 
relatively small n). This figure gives an illustrative sum- 
mary of the physics of the Kondo-lattice model. 

The different behaviors shown in the last three figures 
are worth analyzing in more detail. At high tempera- 
tures, the magnetic susceptibility is decreasing in 1/T as 
normally expected. Then, below the highest character- 
istic temperature (Tcor or Tk), there is a rapid decrease 
of the product xT, indicating a reduction of the effective 
magnetic moment due to the occurence of short-range 
antiferromagnetic correlations, which tend to reduce the 
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response of the localized moments to an applied magnetic 
field, and an even faster reduction due to the Kondo ef- 
fect. At lower temperatures, within the Kondo regime, 
we observe a flattening of the effective magnetic moment. 
This region of nearly constant magnetic moment below 
Tk can be related to the "exhaustion" problem, because 
there are not enough electrons to screen the localized 
spins. At still lower temperatures, a coherent regime 
sets in, and this residual moment becomes completely 
screened. This behavior is not observed in the half-filled 
case, where the effective magnetic moment tends very 
rapidly to zero due to the existence of a gap at the Fermi 
level. 



V. CONCLUSIONS 

We have studied the different behaviors of the intra-site 
and inter-site correlations, and the characteristic tem- 
peratures as a function of the intrasite Kondo exchange 
coupling, of the intcrsite exchange coupling, and of the 
number of conduction electrons relative to the number 
of / electrons, which remains fixed to 1. We use here 
a mean-field approximation with two order parameters 
in each case, the mean-field Kondo correlator and the 
short-range magnetic correlation. 

We, have established here for n < 1, as in previous 
worUD for n = 1, that the dependence of the Kondo tem- 
perature with the coupling constant Jk for the lattice 
can be significantly different from the single-impurity 
case. This result can account for th.e pressure depen- 
dence of Tk observed in|-peRh2Si2j3 CeRu2Ge2j3'El or 
more recently Ce2Rh3Ge5J3 On the other hand, depend- 



ing on the relative values of Jh and Jk, as well as 
on the band-filling, the lattice Kondo temperature can 
closely follow the single-impurity one, as observed in 
many cerium and all ytterbium compounds. Further ex- 
periments are needed to better understand the condi- 
tions yielding a Kondo temperature for the lattice much 
different than the single-impurity one. This issue has 
also been addressed by different theoretical approaches 
to VintVi fViP Knyirln lai-1-ir-p and the Periodic Anderson 

Another interesting result concerns the derivation of 
a correlation temperature below which short-range mag- 
netic correlations appear, in good agreement with neu- 
tron scattering experiments in cerium compounds. These 
correlations can coexist with the Kondo effect, but even- 
tually dominate, and suppress the Kondo regime for suf- 
ficiently high values of the intersite exchange interaction 
or sufficiently low band fillings. 

Our present calculation addresses again the difficult is- 
sue of the nature of the ground state and screening in the 
Kondo-lattice problem. We have shown here that, as the 
number of conduction electrons is reduced, exhaustion 
may be compensated by formation of intersite singlets 
of localized spins. Exact calculations for small clusters 
would be interesting in that a comparison between re- 
sults with a number of conduction electrons equal to, or 
much smaller than the number of 4/ electrons localized 
on the differetit sites of the clusters could shed new light 
on this issue.Ej Finally, it is interesting to notice that tak- 
ing into account lattice effects is essential for describing 
the properties of cerium or other anomalous rare-earth 
compounds at low temperaturpspas shown, for example, 
in photoemission experiments.r'rr^ 
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